Reentrant layering in rare gas adsorption: preroughening or premelting? 
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The reentrant layering transitions found in rare gas adsorption on solid substrates have conflict- 
ually been explained either in terms of preroughening (PR), or of top layer melting-solidification 
phenomena. We obtain adsorption isotherms of Lennard- Jones particles on an attractive substrate 
by off lattice Grand Canonical Monte Carlo (GCMC) simulation, and reproduce reentrant layering. 
Microscopic analysis, including layer-by-layer occupancies, surface diffusion and pair correlations, 
confirms the switch of the top surface layer from solid to quasi-liquid across the transition temper- 
ature. At the same time, layer occupancy is found at each jump to switch from close to full to close 
to half, indicating a disordered flat (DOF) surface and establishing preroughening as the underly- 
ing mechanism. Our results suggest that top layer melting is essential in triggering preroughening, 
which thus represents the threshold transition to surface melting in rare gas solids. 



Rare gas solid surfaces and films provide an important 
testing ground for a variety of surface phase transitions. 
Surface melting 0, roughening and more recently 
preroughening (PR) |^ have been identified or at least 
claimed at the free rare gas solid-vapor interface. Layer- 
ing transitions of thin rare gas films on smooth substrates 
have given rise to a wide literature Q . The discovery of 
reentrant layering (RL) - the unexpected disappearance 
and subsequent reappearance (well below the roughening 
temperature) of layering steps in adsorption isotherms on 
smooth substrates [^,^ - has led to a debate 0| . One pos- 
sible explanation is PR, a phase transition which takes a 
surface from a low-temperature "ordered flat" state, with 
essentially full surface coverage (T < Tp^), to a high 
temperature "disordered fiat" (DOF) state, with half 
coverage, and a network of meandering steps (T > Tpji). 
Layering would disappear at PR, but re-enter in the DOF 
state l^l). The competing explanation is based on the 
possibility of a melting-solidification- melting sequence in 
the top surface layer, similar to that seen for increas- 
ing temperature in canonical molecular dynamics simula- 
tions 1^. In this picture, RL would result directly from a 
layer-promotion-driven melting of the top surface layer, 
and the subsequent advance of a solid-liquid interface 
0. Both approaches appear to capture some important 
physics, but both also have problems. The non-atomistic 
statistical mechanics lattice models provide, in presence 
of an attractive substrate potential, an overall adsorption 
phase diagram with zig-zag lines of heat capacity peaks 
(whose behavior has been called "zippering" |jl^) which 
are centered at Tpn and strikingly resemble experimental 
observations Because they contain PR, the models 
can naturally explain why the coverage jump across RL 
should be about half a monolayer, as seen in ellipsome- 
try Q and in X-ray measurements [Tl| . However, they 
fail to account for continuous atom dynamics, in partic- 
ular melting, and it remains unclear how bad the total 



neglect of these aspect might be at these relatively high 
temperatures. In Ar (111), RL takes place near Q%K, 
not too far from melting at r,„ = SAK . By contrast, the 
atomistic canonical simulation approach does not suffer 
from that problem, and can describe quite well all the 
surface degrees of freedom, including thermal evolution 
of each surface layer from solid to liquid. It finds, realis- 
tically, that top-layer surface melting seems to be setting 
on precisely near the RL temperature. However, it does 
not explain the half layer coverage jump across RL. A 
crucial underlying difficulty of this approach lies in the 
fixed particle number - a difficulty which the lattice mod- 
els, being naturally grand canonical, do not encounter. 
In this situation Grand Canonical Monte Carlo (GCMC) 
atomistic simulation should be the method of choice, ap- 
plied since long ago ||l2| to describe adsorption, albeit of 
a single monolayer. Recently, we demonstrated how a 
free rare gas (111) surface can be realistically simulated 
by GCMC with a Lennard-Jones potential, and found 
indications that PR is indeed incipient at 0.8 Tm [ p^ . 
That work however remained incomplete, because a full 
equilibrium stabilization of the grand canonical surface 
proved to be increasingly hard with increasing temper- 
ature, and failed above 0.8Tm, where a value of that 
would cause neither decrease nor increase of the total 
particle number could no longer be found. 

In this Letter we present results of a fully equilibrated 
realistic GCMC simulation of multilayer rare gas adsorp- 
tion on a flat attractive substrate. In this case, the sub- 
strate potential naturally provides the necessary stabi- 
lization for the system. We obtain realistic adsorption 
isotherms, whose main features compare directly with 
experiment. Reentrant layering is recovered, and layer 
occupancies confirm its association with a DOF surface 
and thus with PR. At the same time however, surface 
diffusion and pair correlations show that while the vir- 
tually full monolayer below TpR is solid, with only a gas 
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of adatoms and vacancies, the half-full monolayer found 
above TpR is made of a 2D liquid islands (even if in a 
strong periodic potential). A new picture emerges, where 
the fractional monolayer melting, besides opening the 
way to surface melting, is also a key element favoring 
the preroughening of these surfaces. 
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FIG. 1. Calculated adsorption isotherms for ti 
(0.53e) and t2 = 0.86 (0.61e), respectively below and above 
reentrant layering {t — 0.83). Arrows indicate that at ti the 
layering steps lead to roughly integer coverage, while at t2 
they lead to roughly half-integer coverage. The abscissa for 
isotherm t2 is shifted by 0.4. Inset: amplitude of the plateau 
breadths between two jumps in the t2 isotherm. Circles: this 
simulation. Squares: extracted from Ar/graphite data 
Full line: model fit to experimental breadths, as described in 
the text, with c/k = 1200. Dashed line: prediction of the 
same model, for a 10-fold enhanced substrate attraction 7. 

We simulate adsorption by classical GCMC, imple- 
menting small displacement moves (m), creations (c), and 
destructions (d) with relative probabilities a*-™-* — l — 2a 
and a^'^' = a^''^ — a. Small moves apply to all parti- 
cles, whereas creation/destruction is restricted to a fixed 
surface region, about four layers wide, since their ac- 
ceptance in the fourth layer of this region is already 
negligible on the entire MC run. In standard bulk 
GCMC the fastest convergence to the Markov chain is 
for a = 1/3 ijl^. For our surface geometry and our 
potential, the optimal value of a is found to be small, 
of order 10^"^ (the precise value depending on the outer 
layer population relative to the total), as needed to al- 
low for a more effective equilibration after each cre- 
ation/destruction move. Creation and destruction ac- 
ceptance probabilities were checked explicitly to satisfy 
the detailed balance. We simulated adsorption of atoms 
interacting via the (12,6) Lennard-Jones potential trun- 
cated at 2.5(7. The bulk fee triple point temperature Tm 
of this model is ^ 0.7e [|l5| (note that pressure depen- 
dence is negligible, i.e. Pm/T„i{dT,n/dP) ~ 2 x 10"'' for 
Ar), and we will from now on switch notation to a re- 



duced temperature t = T /Tm. The substrate was taken 
to be flat and unstructured. Periodic boundary condi- 
tions were assumed along the x and y directions, with a 
reflecting wall along z, placed way above the surface. In- 
teractions between atoms and substrate were also of the 
Lennard-Jones form, giving rise to a laterally invariant 
(3,9) potential V{z) = A{B / z'^ ~C / z^), with A = 407r/3, 
B = 1/15 and C = 1/2, the latter ~ 10 times larger than 
the true Ar/graphite value, so as to avoid the stabiliza- 
tion problems encountered previously with the free solid- 
vapor interface |l3| . The {xy) simulation box size was of 
22 X 23 CT units and a full fee layer contained Ni — 480 
atoms. We focused on two temperatures, ti = 0.75 and 
t2 — 0.86, (respectively below and above the RL temper- 
ature t ~ 0.83), where we obtained full and converged 
adsorption isotherms. For each temperature we increased 
the chemical potential ^ (i.e., increased the pressure of 
the fictitious perfect gas in contact with the system) by 
intervals of ~ 0.02e and waited for stabilization of both 
total energy and particle number. Generally half a mil- 
lion Monte-Carlo (MC) moves/particle were sufficient to 
reach equilibrium. Then 30 to 50 uncorrelated configura- 
tions were generated from a subsequent half million MC 
moves and analyzed. 

Fig. 1 shows the calculated adsorption isotherms - the 
number of adsorbed layers versus (/zo — ^)~^/'^ - no being 
the saturation chemical potential (where a bulk quantity 
of matter would condense). At the lower temperature 
ti we find clear layering steps between consecutive inte- 
ger layers numbers. Analysis of layer occupancies shows 
that after each coverage jump the first layer is nearly full, 
with ~ 15 — 20% of vacancies, and only few adatoms. In 
the subsequent plateau the adatom population gradually 
increases to ~ 15 — 20% and vacancies in the first layer 
are filled, until the next jump suddenly occurs, and so 
on. Between ti and t2 we generally observed that, as 
in experiments, the layering steps tended to disappear; 
however here it became very difficult to obtain a stable 
surface and thus well defined adsorption isotherms. At 
the higher temperature t2 we did recover stability, and 
we found that layering was again present, but with two 
important qualitative differences with the low temper- 
ature isotherm: coverage was shifted by half a mono- 
layer, and plateaus were broader. Adsorption began at a 
half-full layer here, and it progressed continuously, lead- 
ing to a broader plateau, until the next jump to another 
half integer coverage. We plot in Fig. 1 the t2 isotherm 
up to eight adsorbed layers, the maximum thickness be- 
fore encountering again stabilization problems. The large 
plateau breadths are clearly due to our strong substrate 
potential. The film grand potential can be crudely mod- 
eled as a periodic part, say A:cos(27rn), plus an effective 
interface repulsion, c/(2(n — Uq)^) [Q, plus a growth 
term /xn (n is the total number of layers). The plateau 
breadth is thus predicted to decrease asymptotically in 
the form An ~ l/(l-|-7~^(n— no)^), where 7 = 3c/ (47r^A:) 
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measures the strength of the substrate. As Fig. 1 (in- 
set) shows, this law fits well the experimental data, with 
c/k — f200. It also agrees fairly well with our actual 
GCMC plateau widths, once 7 is increased by the cor- 
rect factor 10. We also note from Fig. 1 the relatively 
large compressibility k^^ of the half-coverage state with 
respect to the low-temperature state. 
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FIG. 2. Occupancies of the different layers versus the 
chemical potential calculated at t2 — 0.86 (0.61e). White dia- 
monds indicate a liquid-like layer; black diamonds a solid-like 
layer. 

We conclude that our simulation reproduces the ba- 
sic RL phenomenon, making it possible to probe deeply 
into its nature. For a better understanding of the lay- 
ering reentrance, we plot in Fig. 2 the occupancies, at 
t2 = 0.86, of the different layers for increasing chem- 
ical potential. The jumps leading to fractional cover- 
age states are clearly observable. Following each jump 
(A, layer #5), the coverage increases continuously by a 
fraction of monolayer, enriching the adatom population, 
as well as first and second layers, until at (B) the sur- 
face (layer #6) is ready for the next jump, leading to 
(C) where, following the jump, former adatoms (layer 
#6) increase in density to form a new half layer, and a 
new adatom layer (#7) is started. We found no trace of 
the non-monotonic occupancies reported in earlier canon- 
ical studies Q . The top layer occupancy extrapolates to 
about 50 % for large adsorbate thickness, strongly sup- 
porting the identification with a DOF state: an ordinary 
2D liquid should display a much higher average lateral 
density. The occupancies of the three outermost layers 
(0.1, 0.5, 0.8) for what we thus suppose to describe a 
realistic DOF state differ somewhat from the simplistic 
ones expected from lattice models, namely (0.0, 0.5, 1.0). 
The finding of a DOF surface at against an ordinary 
fiat surface [occupancies (0.15, 0.85, 1.0)] at ti indicates 
that PR of the free rare gas solid surface must take place 
in between. This conclusion is also supported by the ev- 
idence of DOF phase separation taking place at t ~ 0.83 
independently obtained by canonical simulations of the 
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FIG. 3. Pair correlation functions of the states described 
in the text, a) Arrows indicate shell features of a solid layer; 
b),c). Absence of shell features indicates liquid layers. 

One might thus be led to think that apart from details, 
the physics is just that dictated by simple SOS models 
However, a closer look at our MC configurations 
reveals that the situation is different, and richer. Fol- 
lowing ^ we studied the lateral positional ordering and 
diffusion coefficients of different layers at the two temper- 
atures by examining pair correlation functions, in partic- 
ular at t2 = 0.86. For this purpose we carried out two 
separate canonical molecular dynamics simulations (the 
diffusion coefficient is ill-defined in a grand canonical sim- 
ulation), one with an integer layer number, and another 
with half-integer (no substrate) . They were meant to ap- 
proximate free stable grand canonical surfaces below and 
above tjiL = 0.83, and thus chosen with the same cover- 
ages ~ 0.5 and ^ 1 of the grand canonical states A and 
B described earlier. Fig. 3 shows a selection of lateral 
pair correlation functions g{r) calculated at t2- Presence 
of shell-related peaks/shoulders indicates a solid layer, 
their absence a liquid layer . We see that the top layer 
is always liquid, but that it solidifies right after being cov- 
ered by the next half layer. Consider for instance state 
A in Fig. 2. The upper layer (#6) has 10% of adatoms 
(a 2D gas), the lower layer (#4) has 20% of vacancies 
and is solid, but the the middle half filled layer (#5) is 
liquid. As coverage increases, layer #5 gets denser, but 
remains liquid until jump B (see Fig. 3). After that, at 
C, the former adatoms condense into another fiuid half 
layer #6, while at the same time layer #5 solidifies, lead- 
ing to a surface identical to the starting one except for 
one extra layer. This picture is close to that suggested 
by heat capacity studies . It is also similar to that de- 
scribed by canonical simulations 1^ , differing however in 
two crucial respects, namely (i) the lack of a solid-fluid- 
solid evolution for any layer and, more importantly, (ii) 
the half occupancy of the fluid layer. The latter is the 
hallmark of the DOF state, which here therefore emerges 
as the likeliest explanation for RL. 
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In order to further elucidate the connection between 
surface melting and PR, we examined the lateral diffu- 
sion coefficient layer by layer. Mean square displacements 
were averaged for all particles spending time within 
three vertical windows corresponding to adatom layer, 
first layer (surface) and second layer. Not surprisingly, 
adatoms are very diffusive (gas-like), while buried layers 
are solid, and poorly diffusive. The top layer diffusivity 
was always sizable, but larger by about a factor two in the 
half-covered case, where it is similar to the surface mass 
transport coefficient near Tm (Fig- 4). This confirms 
that a height jump by about half a layer across PR also 
takes the top layer from solid to liquid, in agreement with 
the GCMC analysis. Thus sudden formation of the liq- 
uid half layer at PR represents the threshold for the first 
appearance of the liquid, which will subsequently extend 
to lower layers and grow critically to a thicker liquid film 
as temperature is further raised to approach T„i. 
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FIG. 4. Lateral diffusion coefficient obtained for a free 
Ar(lll) surface with half (black and white squares) and full 
(black squares) coverages at different temperatures. The RL 
onset temperature region Tpr is dotted. At half coverage, the 
top layer is always liquid; at full coverage, the top layer is still 
solid at TpR. The second layer is always weakly diffusive, and 
adatoms (not shown) diffuse as a gas in every case. Solid lines 
guide the eye from full coverage below TpR to half coverage 
above TpR. Note the diffusion coefficient jump. Surface mass 
transport coefficient near Tm (from [Q) is dash-dotted. 

Summarizing, our results can explain the experimental 
evidence of RL occurring in the adsorption of rare gas on 
a solid substrate. Layer by layer occupancies and direct 
insight on the surface processes that are not directly ac- 
cessible from experiments, confirm the interpretation of 
the reentrant layering transition in terms of prerough- 
ening. The DOF state consists of a half monolayer of 
barely percolating 2D liquid islands, floating on top of a 
solid substrate. We found a coincidence of the onset of 
premelting in the top layer with a PR transition, where 
coverage jumps from full to partial. These two surface 
phenomena, apparently very different, appear here to be 
intimately connected. A lattice model addressing this 
connection has been published separately . 
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